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Abstract 

The probabilistic Glauber formula for nuclear absorption used in the lit- 
erature is only valid at low energies and x+ ~ 0. Due to energy conservation, 
cr a bs is replaced by an effective cross-section a a b s + x\ (ct C c-n — &abs) which 
increases with x+ and tends to the total cc — N cross-section a C c-N- Experi- 
mental data can be described with a a b s ~ 4 5 mb and <j C c-n ~ 15 20 mb. 
At high energies, due to the increase of the coherence length, this formula 
changes. The main change is the replacement of a a b s by <j C c-n ~ for all values 
of x + , as s — > oo. Thus, if cr C c-N > Cafes the J/ifi suppression due to nuclear 
interaction will increase with energy. 
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LOW ENERGY PROBABILISTIC FORMULA. Let us consider a proton- 
nucleus collision. The probabilistic Glauber formula for the J /if} survival is well 
known. A cc pre-resonant state is produced inside the nucleus at some value Z of the 
longitudinal coordinate, and interacts with the nuclei on its path. If the interaction is 
inelastic, the cc pair can transform into another pair which has vanishing projection 
into J/tf>. The corresponding cross-section is called absorptive cross-section a abs . 
We can write a abs = (1 — £)a C c-N, where a c5 -N is the total cc — N cross-section 
and ea c5 -N is the contribution to <j C c-n of those intermediate states which have a 
non-zero projection into J /if). In the Glauber model one has 



r+oo rca 

a pA( b ) = a pN A dZ p A (b,Z)exp -a abs A dZ'p A (b,Z') 

J —oo L J z 

= (^piv/^) [1 - exp (-a abs A T A (b))} . 
Note that for open charm production a abs = and a pA = Aa pN . 



ASYMPTOTIC FIELD THEORY FORMULA. At high energies the coher- 
ence length increases and the projectile interacts with the nucleus as a whole 0. 
Thus, it is no longer possible to consider its collisions as "successive" and expression 
(HD is changed. As we shall see below this change consists in the replacement [p]] 



[l/a abs ) [1 - exp (-a a bs A T A (b))] =>- A T A (b) exp 



\oA T A (b) 



(2) 



where a = <j C c-n- The change is twofold. There is a change in the form of the 
expression and a a bs has been replaced by <j C c-n- If o'abs ~ o'cc-n the change @ is 
not important. Indeed the two expressions coincide at the first and second order in 
a and, since a abs is not large, the result will not be significantly changed. However, 
if (t c5 -n ^ &abs-, there will be a very important increase of the so-called normal J/if> 
suppression (which can no longer be called nuclear absorption since it is due to the 
total cross-section cc — N, rather than to its absorptive part). This can be the case if 
the cc pair is produced in a color state and is accompanied by light quarks in order to 
have a colorless system. In this case, <j C c-n can be a typical hadronic cross-section, 
due to the presence of the light quarks ||. Furthermore, the probability for this cc 
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system to keep its projection into J/ip after interacting inelastically has to be large, 
i.e. £ ~ 1. 

The derivation of the asymptotic formula is straightforward. Let us denote by 
a the cross-section for the interaction of the light quarks in the projectile and 
a = cr C c-7v that of the produced cc system (we assume, for simplicity, that the 
corresponding amplitudes are purely imaginary). Let us single out one of the light 
quark interactions in which the cc pair is produced (for instance via gluon-gluon 
fusion). We have |L] 

°m( 6 ) = <$v + cO -<r A {J]/a (3) 

where 



a A (a) = 2 jl - exp a A T A (6)] } (4) 



is the Glauber total cross-section in terms of the cross-section a for the scattering 
on a single nucleon. In eq. (|3|) we have substracted the term with no light quark 
interaction (since we assume that at least one such interaction is needed to produce 
the cc pair). The division by a is due to the fact that one such collision is precisely 
the one producing the cc-pair - first factor of (]3|). From (|3|) and M) we get 



<4a$) = WpnIo) ^4(<x)exp 



\*A T A (b) 



(5) 



We see that the terms containing a and a factorize. The former produce nuclear 
effects (shadowing) in the nucleus vertex function (which is no longer proportional 
to A). The terms containing a correspond to the rescattering of the cc system and 
have the form given by (E3). 

As stated above, the main change between the probabilistic formula and the 
asymptotic one consists in the replacement of a a bs by u C c-n- This change is entirely 
due to coherence effects. At low energy, some contributions vanish due to the fact 
that there is a non-vanishing minimum momentum transfer (t m i n ^ 0) - and this 
contribution is suppressed by the nuclear form factor. It turns out that the other 
contributions cancel with each other, giving rise to an A 1 behaviour typical of open 
charm production. The only case in which this cancellation is avoided is when the cc 
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pair is "destroyed", (i.e. converted into a cc system with vanishing projection on the 
J/ip)- For this reason only the absorptive part of the total cc — N cross-section plays 
a role at low energy. At asymptotic energy t min = 0. Thus, the above cancellation 
does not occur and the total cc — N cross-section contributes to the J/tp suppression. 

The expression of the shadowing corrections in eq. (§) is not realistic since the 
Glauber expression has been used. It is, indeed, well known that shadowing cor- 
rections have to be described by triple Pomeron diagrams. In view of that, in the 
following we shall only consider in eq. (|3]) the terms linear in a, and will use a stan- 
dard formulation of the shadowing corrections (which correspond to higher order 
terms in a). It is, nevertheless, instructive to see from the complete expression (|3]) 
that the shadowing corrections vanish in the low energy limit and are factorizable 
in the high energy one. In the following we will assume that factorization also holds 
at intermediate energies - and will denote the corresponding cross-section a^fi a ow _ 

INTERPOLATING FORMULA. Using the AGK cutting rules || it is possible 
to obtain the exact formula at finite energies. We have U instead of (§) 

A / A\ n 

oUW = {°X adow /°) E L E T n ] °$ ( 6 ) 

n=l \ n J j=l 

where 



*<P(b) = o 



a 



3-1 



a 



(-Ifo-D 



-(l-e)a 



3-1 



j[-(l-e)a 



n-3 



3-2 



(j - 1) [-(1 - e)a\ 



n-j+l 



(7) 



and 



T^\b) =n\ dZj dZ 2 --- dZ n cos(A(Z 1 - Z 3 ) ]J p A (b, Z t ) . (8) 
Here A is the inverse of the coherence length 



1 M - 

A = j = m p -££ (9) 

l c S X\ 
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with xp = x\ — x 2 and x\x 2 s = M cg . 

In (H) the longitudinal coordinates of the interacting nucleons have been ordered 
as Z\ < Z 2 ■ ■ ■ < Z n . The index j in (0) and (^) denotes the first interaction 
(i.e. the one with the smallest value of Z) which has been cut in such a way that 
tmin = for the n — j ones with Z > Zj. For details see ref. |J. The damping 
factor cos(A(Zi — Zj)) depends on the difference Z\ — Zj. The only case in which 
tmin = is when j = 1, i.e. when the first interaction is cut. In all other cases 
train 7^ 0. Therefore, the term with j = 1 is the only one that survives in the limit 
s — > 0(A — > oo). It is easy to verify that one recovers in this way the probabilistic 
formula (|]) - which depends only on a a b s - Likewise, it can be verified^] that, in the 
limit s — > oo(A —>■ 0) we recover the asymptotic expression @ - which depends 
only on <t C c-n- 

We see from @ that at fixed s and x 2 , A increases with x\. Therefore, if a C c-N > 
&abs, J/i> suppression will increase with x\ - a feature observed experimentally. In 
ref. an attempt was made to describe the data in this framework. However, as 
recognized in ||, there is a caveat. Since A is a function of x 2 , one obtains a scaling 
in x 2 , whereas the data indicate rather a scaling in xi (or xp). A way out is to take 
into account the modifications of the AGK rules resulting from energy conservation 
H. Indeed, by energy conservation, the x distribution of J/ijj gets softer when the 
number K of inelastic collisions increases. This has been taken into account in |3j 
introducing a factor 

F k (x 1 ) = (l-xJ) k (10) 

in the contribution to a^ A corresponding to k inelastic collisions of the cc system. 

In ref. this has been done in the probabilistic Glauber model. In order to do it 

in the general framework described above, one has to decompose eq. (|7|), into a 

sum of terms a J n k corresponding to a fixed number k of inelastic collisions. It turns 

out that this amounts to replacing [— (1 — e)a] n ~ j by J2 [ n k J )(^) k (—^') n ~ : '~ k Fk{xi). 

With the form of F k in eq. (|10D the summation in k can be performed analytically. 

2 For A = 0, Tn = T\- Moreover, if one considers the two terms in the last equality of (0), 
there is a cancellation between the first term of and the second term of a„ ■ The only term 
left is thus the first term of o^f* . 
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The final result is 



*s°(&) = (-) 3 ' 1 j i-°e ff r j - (-§ )'~ 2 (j - 1) [-<x e// r j+1 (11) 

where 

cr e // = 5(1 — 5) + 5e xj . (12) 

Thus, as a consequence of the modification of the AGK rules due to energy 
conservation, we have obtained the same expression (|7]) with the absorptive cross- 
section a a b s = (1 — e)5 replaced by cr e ff- Therefore, the asymptotic limit, eq. 
which does not depend on a a b s is not changed. 

As discussed above, in the low energy limit, A — > 00, only the term j = 1 
survives. We obtain in this case 

= t 1 - ex P H*// A T *( b ))] ■ (I 3 ) 

We see that for x\ — > we recover the probabilistic expression ([]]) with a e ff ~ cr a 6s- 
However, for 07 — > 1 we have cr e // ~ cr C c-v- While coherence effects produce this 
change from a a b s to cr c5 _Ar as s — > 00, the introduction of the factors Ffc(xi), eq. 
fllOD , leads to the same effect at low s as X\ increases. We recover in this way scaling 
in xi at low energies. 

As shown in ref. [|], eqs. (|P2])-(|l3D give a good description of experimental data 
with a abs = 5 mb, a C c-N = 20 mb (e = 0.75) and 7 = 2. It was also shown in 0] 
that 7 increases with the mass of the heavy-quark pair and thus, the effect of energy 
conservation via the factors Fk concentrates to larger values of x\ when this mass 
increases. 

CONCLUSIONS AND OUTLOOK. We have shown that the probabilistic Glau- 
ber formula for the survival of J/ip in a pA collision changes with energy due to 
coherence effects. The main change is the replacement of a a bs by a c5 -v - the total 
cc — N cross-section. The same change occurs at low energy when x\ increases from 
to 1, due to energy conservation. If a abs ~ ct c5 -n the discussion in this paper is of 
little phenomenological relevance. Nuclear absorption will be similar at RHIC and 
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SPS energies. Moreover, the increase of the J/ip suppression with xp cannot be due 
to the mechanism described here. If, on the contrary, a c5 -N is substantially larger 
than a a b s , the increase with xp is well reproduced and a huge increase of the J/ip 
suppression in pA collisions at RHIC and LHC is predicted. 



Numerical calculations using eq. (|TT|) ar e in progress |J. They should provide a 
more precise determination of the value of <j C c-n from available data. They should 
also allow to determine quantitatively how the xf scaling is converted into a scaling 
in X2 as energy increases - first at low xf and then extending to larger values of xf 
as the energy increases. 
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